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Abstract 

We apply the recently established connection between nonlinear fluid dynamics and 
AdS gravity to the case of the dyonic black brane in AdS4. This yields the equations 
of fluid dynamics for a 2 + 1 dimensional charged fluid in a background magnetic field. 
We construct the gravity solution to second order in the derivative expansion. From this 
we find the fluid dynamical stress tensor and charge current to second and third order in 
derivatives respectively, along with values for the associated transport coefficients. 
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1. Introduction 



Work over the past year has established an elegant connection between long wave- 
length disturbances of asymptotically AdS black brane solutions and relativistic fluid dy- 
namics. The relativistic Navier-Stokes equations emerge from Einstein's equations applied 
to metrics with slowly varying thermodynamic parameters In the context of the 

AdS/CFT correspondence, this construction also provides an efficient means of computing 
transport coefficients in the dual CFT. The details of this connection have been worked 
out for pure gravity in arbitrary dimensions , Einstein-Maxwell theory in 4 + 1 di- 

mensions H[[|,|8| , and gravity coupled to a scalar field in 4 + 1 dimensions §J . There is also 
an extensive earlier literature on the study of transport coefficients within the AdS/CFT 



correspondence; see |I(J for a review and further references. 

In this paper we consider dyonic black brane solutions in asymptotically AdS4 space- 
times. These are solutions of Einstein-Maxwell theory in four spacetime dimensions carry- 
ing nonzero charge densities and boundary magnetic fields. The computation of transport 
coefficients for the dyonic black brane has already received significant attention in light of 
its connection to quantum critical points in 2 + 1 dimensional condensed matter systems 



ll] , [i^ , P^Jl41 , |15| . The approach taken in these works differs from ours, and yields results 
with a different (but overlapping) regime of validity, as we discuss further below. 

Corresponding to the bulk metric and gauge field are a boundary stress tensor and 
current, which obey the conservation lawsi 

V M T 0/x = 0, V M T^ = Be ij Jj , J M = , (1.1) 

where B represents the constant magnetic field on the boundary. The stress tensor is also 
traceless due to scale invariance. Fluid dynamics is a description of the long wavelength 
evolution of conserved charge densities. At vanishing magnetic field, the conserved charges 
are energy, momentum, and electric charge. Also, by "long wavelength" one means long 
compared to the mean free path, which is in turn proportional to the inverse temperature, 
l T ~ =. One then sets up an expansion in terms of the small parameter Zt<9 m , where 
the derivative acts on the conserved charge densities. At nonzero magnetic field the story 
changes. First, there is now a second length scale in the problem, Ib ~ -?=• Second, 
we see from ( |1 . 1| ) that the momentum current J(p*)^ = T lM is no longer conserved. At 
small B there is a separation of the two length scales, It <C Ib, and the appropriate fluid 
dynamical variables depend on the length scale being probed. At intermediate length scales 
It <S I <C Ib, we should treat It9^ as a small parameter but work to all orders in Isd^. 
Also, in this regime momentum is approximately conserved, and so the momentum density 



should be kept as a hydroynamical variable; this is the approach used in HTl| , p!2| , p!3| , p^ , p!5| . At 

3 There is also the dual current, J M = e MQ/3 V a Jp, but this will play no role in our discussion. 
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the largest length scales, / 3> Ib, momentum is no longer an independent hydro dynamical 
variable, but is instead fixed in terms of the energy and electric charge densities. Thus, in 
this regime we can formulate a simpler effective theory by "integrating out" the momentum 
variables. This is the approach we take in this paper. Note that for large magnetic fields, 
such that Ib ~ It, there is no separation of scales and it is clearly inappropriate to include 
the momenta, since they are non-conserved at leading order. Keeping the momentum 
variables in this case is analogous to retaining a very massive particle in a low energy 
effective field theory; it is not incorrect to do so, but it is inefficient. What is being said 
here is not that the momentum densities are being set to zero (they are not), but rather 
that their values are fixed in terms of the charge and energy densities, and so are not 
independent. 



Our goal is therefore to solve the Einstein-Maxwell equations iteratively in an expan- 
sion with respect to derivatives of the energy and electric charge densities, treating both 
Isdfj, and /t<9 m as small expansion parameters. We carry out this task to second order in 
the derivative expansion, following the methods applied to previous examples. This gives 
us the stress tensor to second order. Using the stress tensor conservation equation, we can 
then read off the current to third order in derivatives. This gives us results for a host of new 
transport coefficients beyond what was computed in [ 11 , 12 , [13 , 14] , 15 , including nonlinear 
terms. On the other hand, for the reasons described above, we are restricted to the small 
frequency / wavelength regime, in contrast to |Tl| , |l^ , [T3|JT^ . p!5| where full AC results can be 



found at the level of linear response (at least numerically). Thus, the two approaches are 
complementary. 



Magnetohydro dynamics is of course a highly developed subject, with many important 
applications to astrophysics and plasma physics. Before trying to compare with this lit- 
erature it is important to bear in mind a few points. First, here we are considering the 
fluid dynamics of a conformal field theory, which is quite different than the gas of charged 
particles considered in most applications. Second, most discussion of magnetohydrody- 
namics are in the context of a dynamical electromagnetic field coupled to the charged 
fluid, whereas here we consider a nondynamical external magnetic field. Note, however, 
that that this should be a good model of the dynamical case if the strength of the magnetic 
field is taken to be very large while holding the charge fixed, since then we should be able 
to neglect backreaction from the fluid. 

It should be noted that although we usually discuss the case in which only a magnetic 
field is present and not an electric field, this is not really a restriction since we can always 
apply a Lorentz boost to our solutions to obtain fluid flows with any electric and magnetic 
fields obeying E 2 < B 2 . 
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2. Structure of fluid dynamics 

We begin with a general discussion of the structure of our fluid dynamics, independent 
of the relation to AdS gravity. We consider a general 2 + 1 dimensional conformal field 
theory with a conserved current J M coupled to a constant external magnetic field B. We 
denote the coordinates of our three dimensional manifold x^. 

2.1. Derivative expansion 

There are two conserved quantities in our theory, charge and energy. In local thermal 
equilibrium the corresponding charge and energy densities of our fluid, p and e, are slowly 
varying functions of the coordinates, x M . The current and stress tensor admit an expansion 
in derivatives of e and p, and the equations of fluid dynamics are then given by ( |1 . 1| ) . All 
that remains to specify the theory are the transport coefficients of the various terms in the 
derivative expansions of T^ v and J v . 

By definition, we have 

T 00 = e , J° = p , (2.1) 

with no derivative corrections. ^ To proceed, we need to work out the derivative expansion 
of the remaining components, T° l , T lJ , and J % . 

Our task is simplified by noting that the equations d^T 1 ^ = Be^ Jj imply that the 
spatial components of the current at nth order are related to the stress tensor at order 

77,-1, 

= -ley (doT^-^ j0 + dkT^' 1 ^ . (2.2) 
With this in mind, the remaining equations of fluid dynamics reduce to 

cV + gV i = o , 

d e + <9,T° l = . 
At zeroth order, tracelessness and rotational symmetry imply 



(2.3) 



J>(0)0t _ Q ^ rp(0)ij _ _ £ gij ; — 

At first order we have 

T (1) ° l = (aiSij + a 2 eij)dj€ + (a 3 Sij + a 4 eij)djp , 

rp{l)ij _ Q 



(2.4) 



(2.5) 



4 In discussions of fluid dynamics in which the fluid velocity appears as a hydrodynamical 
variable the charge and energy densities are usually defined in the fluid rest frame. In our case 
there is a preferred frame, defined such that there is a magnetic field but no electric field, and our 
densities are always measured in this frame. 
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where we used ( j2.2| ). The coefficients 04,2,3,4 are functions of e and p but not their deriva- 
tives. 

We next consider terms of second order in derivatives. Here it is convenient to combine 
the two quantities e and p into a two component object 

?o = (e,p). (2-6) 

The second order spatial current is gotten by applying ( |2.2|) to (|2.5|) . Rotational 

invariance implies that T^ 2 -* * can be built out of the following structures: 

d diq a , d q a diq b , t^d^a , d G q a eijdjq h . (2.7) 

Now, in general, if one is working up to nth order, it is permissible to use the zeroth order 
equations of motion to simplify the current and stress tensor at rath order, since the error 
will be of at least order n + 1. The zeroth order equations of motion yield doe = dop = 0. 
Noting that every term in (|2.7| ) involves at least one time derivative acting on either e or 
p, we see that we can set 

T (2)0i = o Thig leaveg T (2)ij ? which hag the structure 

= blld^qaY 1 + b a 2 [e ik d k d 3 q a ] st + bf[d iqa d jqb \ st + bt[e ik d k q a d 3 q h ] st , (2.8) 

where we are using the following notation for the symmetric traceless part of a spatial 
tensor: 

[M tJ ] st = X - (Mij + Mji - 5 %J M kk ) . (2.9) 



Finally, the third order current is obtained by applying ( |2.2| ) to the second order stress 
tensor. 

Our fluid dynamics has thus been specified in terms of the coefficients 04,2,3,4 and 
^1,2,3,4, which are themselves functions of e and p. In the rest of this paper we will see that 
the boundary stress tensor and current obtained from the Einstein-Maxwell equations 
in 4 + 1 dimensions falls into the above framework and yields specific formulas for the 
coefficients. In particular, 01,2,3,4 are given in equation (6.14); b\^ are given in equation 
(6.29); and 63,4 are given in equations (6.34) and (6.35). 

In the above we considered fluid dynamics in a magnetic field but vanishing electric 
field. But by applying a Lorentz boost we can turn on an electric field, subject to the 
condition E 2 < B 2 . As a simple example, consider an infinitesimal Lorentz boost with 
velocity v % : (if = t + v ■ x : x' 1 = x 1 + v l t). The magnetic field is unchanged, and the new 
electric field is Ei = —BeijV^. To keep things simple, let's also work to lowest order in 
both v l and the derivative expansion (i.e. consider v % ~ d^e ~ dip), in which case we can 
write the current after the boost as 
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which displays the Nernst effect as well as the Hall conductivity. We can also write down 
an expression for the energy flux at this order, 

T (1)CH = (axSij + a 2 e ij )d j e + (a 3 <% + a±e l3 )d jP + —e^Ej . (2.11) 

JD 

By performing a finite boost and keeping terms of higher order in the derivative expansion, 
one can of course easily read off the generalized version of these results, if desired. 

We have chosen to formulate our fluid dynamics in terms of the energy and charge 
density, since these have a direct physical significance and appear naturally in the gravita- 
tional description. But we can alternatively work in terms of the temperature and chemical 
potential. The definition of these is subject to some ambiguity in a non-equilibrium con- 
text. One simple choice is to use the same definition between the energy /charge densities 
and the temperature/potential as appears in the equilibrium case, without any derivative 
corrections. The relation between these sets of variables of course depends on the theory 
in question. For the theory corresponding to the dyonic black brane the relations will be 
given in the following section. 

2.2. Normal modes 

We now look for normal mode solutions by solving the linearized equations of motion 
following from (|2.3[). Dropping nonlinear terms, these become 



<9 e + ai V 2 e + a 3 V 2 p = . 

Inserting plane wave fluctuations on top of the constant background, e + § ee lk - x ~ lujt an d 
p + Spe lk ' x ~ lujt , we find 

(-iua 2 -^-k 2 )k 2 5e+(iBuj-ia 4 ujk 2 -^-k 4 )5p = 

2 2 (2.13) 

[ioj + aik 2 )de + a 3 fc 2 5p = . 
Setting the determinant to zero we find two modes 
uj ~ %a\k 2 + ■ ■ ■ , with 5p ~ 

z (ai&£ + a 3 &§) l4 • , s n (2 ' 14) 

u ~ — ^ 2 — - — fc H , with a!Se + a 3 5p^0 . 

2 aiB 



In our gravitational computation we will find a\ < so that the first mode is strictly 
decaying in time. But for the second mode we find that (aib 3 + ; 3 ^ has the same sign as 
3(B 2 — p 2 )a 4 + (B 2 + p 2 )(3B 2 — p 2 ). Since this can have either sign, the /c 4 mode can be 
either decaying or growing in time, depending on the values of p and e. 
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2.3. Comparison with other approaches 

In this section we have developed a formalism for studying magnetohydrodynamics 
which, on its surface at least, represents a significant departure from existing treatments 
in the literature. While our definitions of the transport coefficients are based on following 
consistently the logic of hydro dynamic expansions, we would be remiss not to discuss in 
detail the relation between the different approaches. The purpose of this section is to pro- 
vide that discussion, together with a dictionary between the formalisms where appropriate. 
We conclude this section with a careful treatment of the shear viscosity and its universal 
relationship with the entropy density. 

In most discussions of relativistic conformal hydrodynamics the starting point is a 
Lorentz covariant parametrization of the stress tensor ,i organized in a derivative expansion 
with respect to the fluid velocity as follows: 

(2 15) 

T^ v = v (A^A 1 " 3 (d a Uf3 + d p u a ) - A^cU Q ) , 

with 

A» u = r}^ + . (2.16) 

Here, we have listed all of the transport coefficients present to first order in derivatives of 
the velocity. A particularly good discussion of the relationship between these expressions 
and the dynamics of charged fluids can be found in 0. We note that the charge density 
does not enter explicitly until the next order of the derivative expansion. 

The form ( |2.15| ) can be justified for fluids in Lorentz invariant backgrounds. In par- 
ticular, in thermal equilibrium the fluid can be taken to have any constant velocity, and 
the hydrodynamic expansion then promotes these constants to slowly varying functions. 
The assumption of a Lorentz invariant background clearly does not hold in the presence 
of an external magnetic field; nevertheless, ( |2.15| ) is taken as the starting point in other 
treatments of conformal linear magnetohydrodynamics from AdS/CFT |T3|JT5| , |T6l . As we 
discuss below, it turns out that sensible results can be so obtained provided that one re- 
stricts to the regime of linear fluctuations about thermal equilibrium (which is all that is 
considered in [0,|T5|,|T3] ) , but this approach breaks down in general, as it does not respect 



the principles of hydrodynamic expansions. 

Focussing first on the shortcomings of (|2.15|) , we note that in the presence of a magnetic 



field there is a preferred rest frame, and in thermal equilibrium the fluid will be at rest 
in this frame. Since thermal states are not labeled by an arbitrary constant velocity u M , 



5 In this section we are only going to concern ourselves with the stress energy tensor since the 
conserved current is expressed in terms of the stress energy tensor by equation ( |2.2| ). 
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there is no motivation for introducing a slowly varying velocity field as a hydro dynamical 
variable, which is the basis of ( |2.15| ). On the other hand, one can consider setting up an 
expansion in terms of the deviation of the fluid velocity from its equilibrium value, and 



this is effectively what is done in Jl3 , 15 , 16 



Even more problematic than the above is the fact that (|2.15| ) implies incorrect re- 
lationships between transport coefficients due to an inappropriate imposition of Lorentz 
invariance where there is none. The presence of an external magnetic field explicitly breaks 
Lorentz invariance - there is now a preferred rest frame corresponding to the frame with 
pure B-field and no E-field. However, ( |2.15| ) fails to take this into account and is mani- 
festly Lorentz covariant. Explicitly, the vector -u M = e^ vp F vp breaks the Lorentz symmetry, 
allowing us to separate out the SO (2) vector part u 1 of w M with a projection operator 



'"" '~ 'if + -Lif ( 2 - 17 ) 



\V\ 



With this operator we can separate the transport coefficients according to their SO (2) 
representations. For example, if we look for the structure u % v? in T 1 - 7 and vPvP in T 00 , we 
see that both terms have coefficients fixed by the term eu^u v in (|2.15| ). However, if we use 
( |2.17|) to add additional terms to (|2.15| ) of the form 



AgAj^W 3 , (2.18) 

then the transport coefficient of u l u^ in T ZJ is now independent of the transport coefficient 
of u°u°. 

Interestingly, the imposition of a spurious Lorentz invariance does not cause problems 
at the level of linearized fluctuations about thermal equilibrium. In the preferred frame, 
we write u p = t> M -\-5u^, with 5u° = and keep only terms linear in 5u l . To lowest order in 
this expansion A = A. Since the rest frame of the fluid very nearly corresponds with the 
preferred frame of the B-field, there is no difference between projecting linear fluctuations 
(or derivatives of fluctuations) to either the rest frame or the preferred frame since any 
difference will be second order. At any order in derivatives, SO (2) invariance lifts to lorentz 
invariance for linear fluctuations. This phenomena explains why other authors have not 
run into any trouble using ( |2.15| ): they were considering only linearized fluctuations. In 
order to study manifestly nonlinear properties though, we needed a new formalism. In 
particular, we avoid the problems described above by imposing only the actual SO (2) 
symmetry, and using as hydro dynamical variables only those quantities that can be freely 
specified in equilibrium, namely the energy and charge density. 

With the preceding comments in mind, it is now interesting to compare results ob- 
tained by these different approaches. To bridge the gap between the two formalisms, we 
define the fluid velocity in our formalism by equating the velocity with the Lorentz boost 
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parameter needed to annihilate T . This definition is equivalent to matching ( |2.15|) with 
( |2.5|) and reading off 

1 2 

u l = j—p T ° l = 3^ [(aitij + a 2 e lJ )d 3 t + (a 3 <% + o 4 ey)^p] ■ (2.19) 

This expression is of course corrected at higher orders in the derivative expansion, but 
this version is sufficient for our current purposes. It is worth noting that this velocity 
is manifestly small in the hydrodynamic expansion regardless of the amplitude of any 
fluctuations as it is proportional to the derivatives of our charge and energy densities. 

With our definition of the velocity in hand, we can now use ( |2.15| ) to find a suitable 
definition for the shear viscosity. Since ( |2.15| ) is valid for linear fluctuations, we can use 
(j2.19|) as a dictionary to relate the shear viscosity r/ to our linear transport coefficients. If 
we plug ( |2.19| ) into ( |2 . 1 5| ) and keep only terms linear in our fluctuations then we have 

T (2)v = h [( ai5ik + a 2 e ik )did k e + (a 3 5 lk + a A e lk )d t d k p} st . (2.20) 

Note that a first order contribution in a hydrodynamic expansion with respect to the 
velocities is equivalent here to a second order contribution with respect to the charge 
and energy densities since the velocity itself is now being considered first order in the 
hydrodynamic expansion. As noted above, second order linear contributions to the stress 
tensor from derivatives of the charge density were not included in (|2.15| ), so the second 
half of ( |2.20|) is not predictive of any relationship with our current results. This means 



that we should only interpret the shear viscosity as a factor in the diffusion of energy and 
not charge. Since the only terms contributing to the second order linear energy transport 



coefficients from the standard formalism are those in (|2.20|) , we can read off relationships 
between our transport coefficients and the shear viscosity r\. Comparison with ( |2.8j ) 

yields 

bt = ^, h\ = ^. (2.21) 
1 3 2 3 K ' 

If we skip ahead and insert the values for these transport coefficients that we obtain from 

AdS/CFT in section 6, we see that both of these equations imply that 

r)=^, (2.22) 

where a is related to the entropy density s according to s = ira 2 . Putting these facts 
together yields the universal ratio 

- = ]~ ■ (2-23) 

S 47T 



This result is a slight extension of [16] which demonstrated this result for uncharged fluids 



in the presence of a magnetic field. It should be noted however that while our result takes 
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into account additional degrees of freedom in the form of a variable charge density, we have 
not used any nonlinear effects to demonstrate this relationship. As noted above, the shear 
viscosity is a purely linear property of our system and has no dependence on nonlinear 
properties. 



3. Gravity side 

We now turn to the gravitational description of fluid dynamics. In this section we 
write down the dyonic black brane solution and its thermodynamic properties, and give 
the general prescription for computing the boundary stress tensor and current. 

3.1. Action and dyonic black brane solution 

The bulk Maxwell-Einstein action with a negative cosmo logical constant, which can be 
thought of as a consistent truncation of the theory resulting from an S 7 compactification 
of M-theory, is 

o r 

d A X\ 



2 

^4 



- a R-- a FmnF - m 



(3.1) 



We henceforth choose units with L = 1. The equations of motion are 

Wmn = Rmn + %9mn — ZFMpFff + \gMNFpQF p ® = 
Y N = V M F MN = . 

The dyonic black brane solution in Eddington-Finkelstein coordinates is 

ds 2 = 2dvdr - U(r)dv 2 + r 2 dx i dx i 



(3.2) 



F = A^dr A dv + Bdx 1 A dx 2 



(3.3) 



where U(r) is given by 



^)- 2 + ^-|- 0.4) 



The horizon is located at the largest real root of the equation U (r) = 0. Calling this r = a, 
a is thus given by the largest real root of 



a2 + p 1 +B 1 _2 f = 
er a 



The energy and charge densities (defined below) are given by e and p. 
The Hawking temperature is 



T = 3a _ B 2 +p 2 
•Itt \-<i :i 
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e and p are restricted to values such that T > 0. The chemical potential can be read off 
from the asymptotic value of A v in a gauge such that A v vanishes at the horizon. This 
gives n = £ . 

Some further conventions: Latin indices M, N, ... run over all four spacetime coordi- 
nates, while Greek indices p, v, ... run over the three coordinates (v , x , x 2 ). Since v plays 
the role of time on the boundary, we will sometimes use v = x°. The boundary theory will 
always see a Minkowski metric, ^^dx^dx 1 ' = —(dx ) 2 + dx l dx l . Indices on the boundary 
stress tensor and currents are raised and lowered with this metric. 

3.2. Stress tensor and current 

The action ( |3.1J ) should be supplemented with the boundary terms 

Sbndy = ~\ I d^Xy^O ~ \ I d^Xy/^J . (3.7) 
K 4 JdM K 4 JdM 

Here 7 is the boundary metric and 9 = 7 At "0 Ai „, where 9^ = —\{^J^n u + V^n^) is the 
extrinsic curvature of the boundary, defined in terms of the outward pointing unit normal 
vector n. 

The conformal boundary metric is defined as 7 Mi , = lim^oo ^2-7^. Also, the boundary 
gauge field is defined as linv^^ A M , in a gauge where n M Am = 0. The boundary stress 
tensor and current are then defined as 



6S = \ [ (2 J^A^ + T^8%„) . (3.* 

K A JdM 



Explicitly g7 



= r 2 F »r 

r 5 (3.9) 

T^ v = [9^ v - #7^ - 27 Ml/ ] 

Implicit in ( p.9|) is the large r limit, as well as a projection of T^ v parallel to the boundary 
(since the orthogonal component does not appear in (|3.8|).) 

Electromagnetic gauge invariance implies current conservation, 

V^J^ = . (3.10) 

Invariance under diffeomorphisms generated by vector fields tangent to the boundary yields 
the (non) conservation equation 

V V T^ = F^J V . (3.11) 
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Tracelessness of the stress tensor follows from invariance under diffeomorphisms shifting 
the radial location of the boundary 



lp. v T» v = . 



(3.12) 



In particular, the latter invariance follows from the absence of logarithmic divergences 
in the bulk action, the presence of which would necessitate adding a non-diff invariant 
counterterm |TH . 



Applied to the solution fl3~3] ) we find 



= diag(e, ie, \e) 
J" = (p,0,0) . 



(3.13) 



4. Gravitational derivative expansion 

The dyonic black brane solution above corresponds to a fluid in thermal equilibrium. 
To find a solution corresponding to local thermal equilibrium, we follow the logic in |]J and 
construct a solution which looks locally like the dyonic black brane. We take the expression 
in ( |3.3|) and allow e and p to be functions of the boundary coordinates x". While this 
expression is no longer a solution of the Maxwell-Einstein equations, it is approximately 
a solution in the limit that derivatives of e and p are small. By adding small corrections 
order by order in a derivative expansion we construct a perturbative solution which locally 
approximates the dyonic black hole. 

To implement the derivative expansion, we formally regard e and p as being functions 
of ex", where £ is a formal expansion parameter that is eventually set equal to 1. The 
metric and gauge fields are expanded as 

g = g<® (e, p) + eg^ (e, p) + e 2 g^ (e, p) + 0(e 3 ) 

A = A^(e lP ) + eA^(e lP )+e 2 A^(e lP ) + 0(e 3 ) , * ' 

where g^ and A^ represent the lowest order solution given in ( |3.3| ). The energy and 
charge densities are themselves given by an expansion, 

e = e {0) (ex") + ee w (ex^) + ■■■ , p = p (0) {ex") + ep^ (ex") + ■■■ (4.2) 

As explained in 0, the equations of motion can be solved "tubewise" by working near a 
given x" location, say x" = 0. It is then convenient to set e < - n>0 - ) (0) = p( n>0 \0) = 0. In 
the following, the order of a quantity is defined with respect to its associated power of 
e. Note that this is slightly different than the labelling used in section 2 where we just 
counted the number of derivatives acting on e and p. 
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For the metric and gauge fields we make the gauge choice 

A r = 0, g rr = 0, 9^9^=0. (4.3) 

The zeroth order solution preserves SO (2) rotational symmetry, and this can be used 
to classify the corrections to the metric and gauge fields. The gauge field corrections are 
written 

A (n) = A (n) dv + A in) dx i ^ ^ ^ 

with A^ an SO (2) scalar and A^ 1 an SO (2) vector. For the metric we write 

(ds 2 )H = —dv 2 - 2h^dvdr + r 2 h^dx l dx l + 2j) n) dvdx i + r 2 a^ ) dx i dx J . (4.5) 

In this expansion, JfcW and are SO (2) scalars; j\ n ^ is an SO (2) vector; and affi is an 
SO (2) symmetric traceless tensor. 

We impose the following large r boundary conditions on the n > components: 



(4.6) 



fc(»)~r° /i (rt) ~- 7^~i cr (n) ~l 

These conditions follow from a combination of the asymptotic AdS boundary conditions 
along with the freedom to redefine coordinates as well as the zeroth order solution, as in 

I- 

In addition to these large r boundary conditions, we must also demand that our 
solution be smooth across the horizon at r = a in order to uniquely determine a solution. 



As shown in , this condition is equivalent to demanding the presence of purely ingoing 
modes at the future horizon. 

Although the Maxwell-Einstein equations of motion are invariant under electric- 



magnetic duality, the boundary conditions (|4.6|) are not. In particular, our boundary 



conditions fix a constant magnetic field and vanishing electric field on the boundary. How- 
ever, electric-magnetic duality will turn on a nonzero electric field in general. See |2(J for 
further discussion of the breaking of duality by AdS/CFT boundary conditions. 
We then find the current and stress tensor (a large r limit is implicit): 

r = P 

r = -r 2 ^d r A\ n) 

n 

rpw 

(4.7) 



n 



^ = -e + -r 3 V <j\ 

9 A l - 



(n) 
ij 

n 
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5. Structure of perturbation theory 

The nth order metric coefficients are determined by the components of the Einstein 

(n) 

equations W MN = with M, N / v . These equations can be organized as 

= -±d r (r 4 d r h^) - S$ = 
r 2 (UW rr )^ - Wi n) = d r + d r (d r {r 2 U^)h^) - ^B 2 h^ + 4p^d r A^ - S$ = 

= \rd r (y 2 d r {rj<? >)) + - Be l3 ]d r Af - S \ n) = 

- U.Mi = dr {-l^dMf) - S\f = 

(5.1) 

The source terms appearing above are constructed from the solution at order n — 1, and 
so are assumed to be known. 

Similarly, two components of the Maxwell equations yield 

Y (n)v = }_g r ( r 2Q rA {n) _ 2p (0) h (n)\ _ y{n) = q 

\ ( 1 \ ( 5 - 2 ) 

Y(n)l = ^ f U^drAf + + Be^jfA - Vt ] = 



These equations, together with the boundary conditions (|4.6|) , are sufficient information 
to solve for all the metric and gauge field functions in terms of e and p. The remaining 
Einstein-Maxwell equations then becomes conditions on e and p; in particular, they are 
identified with equations of fluid dynamics (|1.1| ) expanded to a given order in e. 

6. Solving the equations 

6.1. Zeroth order solution 

At zeroth order we use the solution (|3.3|) but with e = e'' -* and p = p(°\ By construc- 
tion, this is a solution to the Einstein-Maxwell equations. The current and stress tensor 
are 

J» = (p,0,0) 

1 1 (6-1) 
= diag(e, -e, -e) . 



2 7 2 



6.2. First order solution 
At first order we write 



e(x^) = e (0) +£a^e (0) (0) 
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The first order sources are built out of d^e^ and d^p^°\ and read 

c(i) _ c(i) _ c(i) _ c(i) _ T/ (i) _ n 

v m = fi^ (6 ' 3) 

i r 4 



All of the functions in (|5.1| )-(|5T2|) appearing without sources can be set to zero, either due 



to the boundary conditions or by absorbing into the lowest order solution. What remains 
is to find d r A^P and j^. 

Integrating Yv 1 )* = gives 



u(o) drA W + l_ [p (0) Si . + /ie .. :i (D + 9_iP_ = _ c (i) ; {6A) 

where Ci is r-independent. Solving for d r A^ and plugging into = gives 

^(0)^.(1) _ umj (D = ±_ {pi o )5i . _ Bei . )(c (D + l djpi0)) (6 5) 

where ' denotes the radial derivative and we used the following identity for U : 

C" = ^ + ^(P 2 + S 2 ). (6.6) 
The general solution of (|6.5|) is 



with integration constants and M • The integral is elementary but lengthy — see 
the appendix. The large r boundary condition fixes = oo. /?j is fixed by demanding 
regularity at the horizon. In particular, for generic pj we find that j!> has a logarithmic 
term, j^' ~ (r — a) ln(r — a), which yields a divergent derivative. We can cancel this term 
by choosing 

w 4_ / gW(gW) \ w 2_ / 2t/-(°)(«(°)) \ (0) 

P * " a<») V a(0)f/"(»)(a(0)) / ) ' a (o) 2 V + a (°>t/»«D(a«V ,P ' 1 ' 

With this choice of integration constants we find that j- has the large r behavior 



_ (,('% - Bc^ | g 1 
3r r 2 



(0 = 3r ^ 3 + O(-) • (6.9) 
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(1) 

Up to first order the current and stress tensor can now be computed as 



d r A\ can now be found from 



r = & 

T* = -±{pV>5 ij -Be ij )pV (6.10) 



[T ij ] st = . 



The remaining Maxwell-Einstein equations are W^j = and y W r = 0. With the 
above expressions for j t . and d r A^ we can verify that these equations are equivalent to 
the order e part of the equations 

8^ = 0, d u T» v = F» V J V . (6.11) 

Note that for the powers of s to match up, the stress tensor and current appearing with 
a derivative are given by their lowest order expressions ( |6.1| ), while the undifferentiated 
current is given by ( |6.1U| ). In particular, these equations imply 

a„p(°) = d v e^ = 

1 (6-12) 

c 1} = -— e-5-e 
2B 

We can now read off the first order transport coefficients by writing 

T m = (a-iSij + a 2 e ij )d j e + (o 3 <5^ + a 4 e ij )d j p , (6.13) 
as in ( |2.5| ). Using the explicit form of U as well as the relation between e and a, we find 

ai = - 



4 \p 2 + B 2 
3 / e \ p 



tt2 4 \p 2 + B 2 J B 

If p \ ( 3ae — p 2 — B 
a 3 = - 



(6.14) 



2 \p 2 + B 2 J V a 2 
1 / B \ ( 3ae — p 2 — B 2 



° 4 2 \p 2 + B 2 J V a 2 
5.5. Second order solution 

e and p are now given by expanding out (|4.2[) to order e 2 , 

e ( x f) = e (o)( ) + £ ^<9 At e ( °)(0) + -e 2 x^<9^e ( °Ho) + £ 2 x M d M e (1) (0) 
p(a;M) = p (o)( ) + ex ^a M p(o)( ) + l e 2 x"x u d fl d u p^(0) + e 2 x fx d fl p^(0) 
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(6.15) 



The second order sources work out to be 



sf> = o 

4 r „ m „ .nw 4 „ r „ .m (iw 



(6.16) 



+ 2[/(°)[5 r ^ 1) a r A;. 1) ]^ + -[d iP ^d r AfY + ^B[d r A^e jk £ 



i^-*2 ^ J ■ 2. 

v z (2) = o . 

We can solve the equations separately in the scalar, vector, and symmetric traceless 
tensor sectors. 

scalar sector : 

We determine from the equation Wr? = 0, which takes the form 

d r (r 4 d r h 2 ) = -2r 2 (d r ylf } ) 2 . (6.17) 

The source on the right hand side is determined from (|6.4| ), (|6.7| ), and (|6.12| ). Of main 
relevance is the fact that d r A^ is a smooth function, including at the horizon, and has 
asymptotic behavior 

(i) i 

d r A^ = -5L_ + . (6.18) 

We can therefore proceed by straightforward integration: 



h^(r) = -2 f dr" (r"d r Af\r"^ " 

J oo ^ J oo 



(6.19) 



Next, 5 r Al 2 ' ) is found from integrating Y( 2 )" = 0. Writing this equation as 
p-<9 r ^— r 2 5 r Al 2 ' ) J + X\ = 0, the solution with the required boundary conditions is 

d r 4 2) (r) = -W ^V 2 Xi(r') . (6.20) 

T J oo 

Finally, is found by integrating r 2 (UW rr )^ — w!j?' = 0. Writing this equation 
in the form d r (— 2 A^ 2 )) + X2 = we can straightforwardly integrate to get 

1 /" r 

jfe( 2 )(r) = -r / dr'X 2 (r') . (6.21) 

2 ioo 
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This completes the solution in the scalar sector. The explicit integrals appearing above 
are complicated and not especially illuminating. Since these functions do not show up in 
the stress tensor and current it is not necessary to compute them explicitly in order to find 
the second order transport coefficients. 

vector sector : 

The vector sector is simple since the corresponding source terms in ( |6.16| ) are zero. 
The solution is obtained as in the first order computation, with the only difference that 
we do not have the analog of the source terms proportional to dip(°'. We thus have 

j?\r) = -I7< >(r)(p<% - fley) Jjr' \^ t) i , (6.22) 

with 

d r A^ is determined by the analog of (|6.4| ), 

U^d r A^ + ![,<% + Bc iL ;jf = -cf . (6.24) 

Pf> and cX> now contribute to the stress tensor and current precisely as in (|6.10| ). Also, 
the stress tensor conservation equation now fixes 

c? = • ( 6 - 25 ) 

just as in ( |6.12| ). Therefore, the vector components of the stress tensor and current at this 
order take the same form as at first order, 

J 1 = — —eadn6 

2B 13 3 (6.26) 

T vl = {a\8ij + a 2 tij)djt + (a 3 5ij + a A eij)djp . 

with the same coefficients ( |6.14| ). In writing this expression for T v% we have used p( 1 '(x fJl ) = 
0. This is consistent, since the equation d^J^ = at order e 2 is simply d v p^ = 0, as 
follows from the fact that the spatial current is divergenceless at this order. 
This completes the solution in the vector sector. 

tensor sector : 

(2) i (2) 

The tensor components are determined by solving — 2^ij^kk = 0- Integrating, 
we obtain 

^ ) «=- 2 £^v)/J/ / ' 5 « ) ^- < 6 - 27 » 
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The lower integration limit on r" is chosen to avoid a singularity at the horizon, while the 
integration limit for r' is chosen to impose the correct asymptotic behavior. 
Of primary interest is the contribution to T %0 given in ([O]). This gives 



[T 



ij].st 



(o) 



dr sfj (r) 



(6.28) 



From this we can read off transport coefficients 61,2,3,4 appearing in 

The simplest contributions are the second derivative terms, bf [didjq a ] and b% [ei k d k djq a ] st , 
which only get a contribution from the source term sjj = —d r [dijj] st . We find 



K 

b% 



b P 2 



1 



o(°) V'(°)(a(°)) 

a(°) 3 £/"(°)(a(°); 
,(o) 



4 (p 2 + B 2 ) 



ap 



2 (p 2 + B 2 ) 



P x 



a 



P 



(6.29) 



a(°) 2 U"(°)(a(°))B 4(p 2 + B 2 )B 

2B _ 1 aB 

a(°) 3 t/"(°)(a(°)) ~ 2 (p 2 + 5 2 ) ' 



The remaining nonlinear transport coefficients are much more tedious to obtain. There 
are two contributions to the stress tensor: [T lJ ] st = [Tl J ] st + [T^ 1 ^ 81 



The contribution from the source term «S',-,- 



(2) 



—d r [dij^] st is easy to obtain, and 



yields the stress tensor 



ij 1 st 



1 d f a(e, p)' 



+ 



+ 



4de \p 2 + B 2 
1 d f a(e, p) : 



t 1 d f a(e,p)p 



2 dp \p 2 + B 2 
1 d f a(e,p)p 



[dipdjp] 



st 



4 dp \p 2 + B 2 
1 d ( a(e, p) 2 p 



2de \p 2 + B 2 

B d 



[e ik d k edje] st + 



[diedjp] 
a(e,p) 



2 9p \p 2 + B 2 



(6.30) 



[e ik d k pdjp\ 



st 



+ 



B d 

+ "7T 7T~ 



a(e,p) 



2 a e Vp 2 + s 2 y 



[e ik d k ed 3 p] 



st 



AB de\p 2 + B 2 
"_L__d_ / a(e,p) 2 p 
4S dp V P 2 + B 2 

Here a(e, p) is obtained by solving 

The rest of the contribution take the form of integrals, which we will not evaluate 
explicitly. The sources are built up out of j t - and d r A± . These were found in the first 
order analysis, and we write them in the form 



d r A % 



(i) 



(pSij - Btij)(j e (r)e jk d k e + j p (r)djp) 
A e {r)tijd 3 t + A p (r)dip 



(6.31) 
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with 



Je{r 



A e (r) 



Mr) 



(r) = ——U(r) I dr' 



Ba 



1 



1 + 



aU"(a) 



U(r') 2 

2U'(a) 



U(r 



l\2 



(6.32) 



1 a*U"{a) 

2BU(r) ~ 4 r 2 U(r) 
1 1 a 4 U"(a) 



je(r) 



rU{r) 4 r 2 U{r) 



In terms of these functions, the contributions to the stress tensor become 



- £ dri (p 2 - 5 2 ) (- 2 M + 2 ^£ + ^ + 2C/A 2 + l^e^e]- 
+ £ dr{(p 2 - 5 2 ) (-2^ + 2^ + IjX) + WA\ + 45 2 ^ + 4^ J[^p]-* 
dr( (p 2 - E 2 ) (_2 Uej ' p+j,Jp) + 4% + + 4EM e A p +4^ 



+ 4pS £ dr (M _ Mi _ ^ + Ml) [e^^ 

9 Joo V r ~ " 2 3Jp + V 2 J [ 1 jf>] ' 

(6.33) 

Adding these to the terms in (|6.32| ), we can read off all of the second order nonlinear 
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transport coefficients, 



_1 d ( a(e,p) 2 

°3 — ~ — 



w 3 



4de Vp 2 + ^ 2 
1 d f a(e, p) 2 \ Id ( a(e, p)p 



rfr<j (p 2 - S 2 ) ( _ 2 M + 2 M + W) + 2 [/A 2 + 45 2 ^f 



4dp\p 2 + B 2 J 2de\p 2 + B 2 J_ 
+ 4pB [ dr 



2 f)3 £ Jp + 



3 2<9pVp 2 + 5 2 



+ / drl (p 2 -B< 



and 



!^ + 2^ + + 2*7Aj + 4B 2 ^+ 4^ 



4 
r 



tee 
°4 



°4 



1 d f a(e, p) 2 p 



45 de V P 2 + # 2 
1 d ( a(e, p) 2 p 



4p5 / dr 



jej' e jeje 1 -/ ./ . ^-eje 
T^eJe + 



r 



r- 



+ 



B d f a(e, p) 



4S dp V P 2 + 5 2 

dr<( (p 2 -B' z ) ( -2 

+ 4S 



2 <9e Vp 2 + 5 2 

Uejp+jUp) , ,jejp 



2 (A J p + A p j e ) 



(6.34) 



+ 4 ^+« +4C/A e A, + 4-^ (6 . 35) 



B d f a(e,p) 



2 <9p VP 2 + 5 2 



+ 4p£ 



dr 



9 A J pJ p ' 9 



The final step in the construction of our second order solution is to check that the 
remaining Einstein-Maxwell equations, W^ v = Y^ r = 0, imply the equations of fluid 
dynamics ( |1 . 1|) at order e 2 . This is indeed the case. 



7. Conclusion 

We have shown how the dynamics of a dyonic black brane can give rise to nonlinear 
magnetohydrodynamic equations of motion. We explicitly solved the Einstein-Maxwell 
equations to second order in the fluid dynamical derivative expansion, and from this we 
found expressions for the stress tensor and current to second and third order in derivatives 
respectively. These expressions yielded gravity predictions for the magnetohydrodynamic 
transport coefficients of charge and energy densities in a specific conformal fluid, coefficients 
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which can in principle be compared with experimental results dealing with 2+1 dimensional 
fluids near a quantum critical point. Our linear transport coefficients extend to higher order 
results already appearing in the literature, and a new set of nonlinear transport coefficients 
is now available as well. We found explicit closed form expressions for the second order 
linear transport coefficients, and integral expressions for the nonlinear ones. 

Our work also helped to demonstrate the broad applicability of the mechanism |l] 
by which the equations of fluid dynamics emerge from gravity at long wavelengths. Com- 
pared to previous examples, one novelty of our setup was that it did not involve considering 
boosted black branes, due to the breaking of Lorentz invariance by the background mag- 
netic field. Instead, we only needed to consider fluctuations of the energy and charge 
density. 

Various generalizations and extensions of this work are possible, such as considering 
higher dimensional spacetimes, curved boundaries, etc. As with other examples, it may be 
interesting to consider in more detail particular solutions, and to see what can be learned 
about gravity from properties of the fluid dynamics, and vice versa. 
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Appendix A. Integrals 

Here we give the results for the integrals appearing in ( |6.7| ). First write U in terms of 
its roots as U = ^(r — r\)(r — r 2 )(r — r 3 )(r — r 4 ). We note that one of the roots is equal 
to a, while the other three can be obtained explicitly as the roots of a cubic equation. 

We will use the identity 

J_ = V — 1 y U "M 1 f Ai) 

V 2 [U'in)} 2 (r - n) 2 ^ [U'{ ri )f (r -rA 



along with 



Y( 1 l + jgfoLi) =0 



( 2 \TP(~ 



1 . U"(rA 1 



+ r t* t" / / =0 . 



[U'(n)] 2 r? [C/'(r,)] 3 r 2 



(A.2) 



The identities in ( |A.2|) following from applying derivatives to ( |A.1| ), evaluated at r = 0. 
We then find 
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f*L V 1 1 V U "W l ( - ) 

J JJ2 ~ [C/ /( r .)]2 (r _ r .) [C //( r .)]3 n ^ ^ 

r dr ^ 1 / 1 ln(r-r,)\ ^ t/"(ri) ln(r 

y rfp - - 2^ [U'(ri)] 2 yr^r-r,) + 7f ) ~ ^ [[/'(r,)] 3 r. 

f dr v 1 / 1 ln(r-r t ) \ _ ^ U"( n ) ln( 
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